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7 Abstract 

8 A new scenario for the interaction of a jet with a background medium with cold 

9 clouds is investigated by means of hydrodynamic turbulence simulations with cool- 

10 ing. The idea is that the cold clouds are overtaken by a radio cocoon and stirred 

11 up by turbulence in this cocoon. The 2D multiphase turbulence simulations contain 

12 all the three gas phases and have a number of interesting properties. The produced 

13 power spectrum is proportional to the inverse square of the wavevector. The Mach 

14 number - density relation may explain the observed velocities in emission line gas 

15 associated with radio galaxies. The model also explains the increased optical visi- 

16 bility within the radio structures, the correlation between emission line and radio 

17 luminosity, and the evolution of the alignment effect with source size. I also propose 
is this to be a useful model for the recently discovered neutral outflows in nearby radio 

19 sources. 
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23 1 Introduction 



Extragalactic jets consist of relativistic, magnetised, radio-emitting plasma. 
They are driven into the warm, usually X-ray-emitting also magnetised gas 
that surrounds galaxies. Colder gas components are observed in many systems 
as well. There is the alignment effect, i.e. optical emission lines and continuum 
co-spatia l with th e area enclosed by the radio structures at redshifts beyond 
z ~ 0.5 ( Longair . 2006| ). Associ ated Lyman a a bsorption sometimes at the 
same redshift as molecular gas ( KrauseL 2005a ) is detected at even higher 
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31 redshifts. Atomic and io nised hydrogen has also been detected in interaction 

32 with local radio sources ( Morganti et al. . l2005al fbh. 



33 Understanding the interaction of these components is not only an interesting 

34 piece of astrophysics in itself, but bears the potential to find out more about 

35 the galactic nucleus that causes the jets, the environment of galaxies, and the 

36 way jets influence the evolution of galaxies. 

37 2 Turbulence in jet cocoons 



38 The alignment effect describes the observation of cold (< 10 6 K) gas associated 

39 with radio structures. A common idea is that the cold gas is present in the 
form of clouds, with low filling factor, which predate the jet phase. Emission i s 
thought to be produced either by shocks in these clouds ( Mellema et allEo02j ). 
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induced by the bow shock, or by photoionisation due to a hidden quasar (there 
is also emission ascribed to newly formed stars). 



44 A problem with this scenario is the sparsity of emission outside the radio 



structures ( Inskip et al. . 2002^ 1 . where the clouds should still be photoionised 



by the quasar. The obvious alternative, cooling the cold gas from the sur- 
rounding warm phase (10 6 K < T < 10 8 K) was so far excluded, since the 
cooling time of that gas, where observed, was much longer than the jet age, 
and also for dynamical reasons. 



We have found a way to circumvent these difficulties. The new idea is that just 

51 a few preexisting cold clouds are overtaken by the jet. When they enter the 

52 radio cocoon, together with surrounding warm gas, they are stirred up by the 

53 violently turbulent plasma in the radio cocoon and form a turbulent system. 

54 Now, turbulence transports energy between the phases, and the cold clouds 

55 radiate it away. Hence, the whole system looses energy, as does the warm 

56 phase. The turbulently enhanced cooling rate leads to cold mass condensation. 

57 The increased cold mass is in turn responsible for an increased visibility of 

58 optically emitting clouds, even when considering photoionisation as ionisation 

59 mechanism. 

eo 2. 1 Simulations 



An idealisation of the proposed scenario is the Kelvin-Helmholtz instability 
with cool clouds in the denser phase. Since radio jet cocoons are usually less 
dense than the surroundings, this may represent the situation at the contact 
surface between radio cocoon and shocked ambient gas with the expanding 
cocoon just about to overtake some cool clouds. We have simulated this setu 
(see Fig. 1 for a sketch) with th e two codes, Nirvana ( Ziegler and Yorkel . ll99 



67 and Flash ([Frvxell et all 120001 ) . The physical model is hydrodynamics (conti- 
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Fig. 1. Sketch of the setup for all simulations. The three phases soon mix due to 
fluid instabilities, and form a turbulent region powered by the initial gravitational 
and kinetic energy. 

es nuity, momentum and energy equations) plus constant gravity and opti cally 

69 thin cooling. We adopted the approach of iBasson and Alexand"er (2003}), i.e. 

70 a cooling function which is non-zero at all temperatures due to molecules, 

71 atomic emission lines and bremsstrahlung in order of increasing temperature. 

72 Common sense suggests that the Flash simulations should be more accurate 

73 due to the higher order of the interpolation scheme and the strict enforce- 

74 ment of energy and momentum conservation. However, practical constraints 

75 (running time and the ability to cope with high density contrasts) made the 

76 Nirvana code the better tool for this particular job. Also important: to have 

77 similar results with different codes in particular with very different diffusivities 

78 gives more confidence in the conclusions. 

79 We start with a 2D-Kelvin-Helmhoitz setup with a density ratio of 10 4 (Nir- 

80 vana) and 10 3 (Flash). The less dense material has a tangential Mach number 
si of 0.8. Here I restrict myself to the discussion of cloud densities of ten and 

82 a hundred times denser than the denser medium, and a control run without 

83 clouds. 



84 2.2 Results 



85 The evolution of the density for the Flash run is shown in Fig 2. The cloud 

se is shocked (the wobbling of the contact surface sends shocks into the upper 

87 medium), compressed, torn apart, and forms little cloudlets that finally drift 

as away from each other. The Nirvana simulations look quite similar with the 

89 exception that they expected, more diffuse and filamentary. 

90 The turbulence we find is fairly isotropic with a power spectrum proportional 

91 to k~ 2 , where k is the wavevector. This is what would be expected for uncor- 
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Fig. 2. Evolution of the logarithmic density for the Flash simulation. Time increases 
(0, 0.4, 0.5, 0.6, 0.8, 2.6 Myr) from top left to bottom right. Two plots are show at 
each time, the lower ones showing a magnification of the region marked with a box 
in the upper ones. The box marks always the same region. 



related shocks ( Elmegreen and Scald . 2004 ). and may be a special property of 
such multi-phase turbulence. 



94 We always find a Mach number - density relation as shown in Fig. 3 (left). For 

95 the low density, non-cooling part, the Mach number is constant. Because the 

96 pressure generally varies very little across the computational domain, this 

97 means that the kinetic energy density is independent of the density. The 

98 denser, cooling part has a Mach number proportional to the square root of 

99 the density. This means that the velocity is roughly constant in this regime. 
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Fig. 3. Left: Typical Mach number - density histogram here for a Nirvana simulation. 
Gas that is dense enough to cool on the simulation time has constant velocity 
(M oc i/p), non-cooling gas has pv 2 = const (M = const). Right: typical gas mass 
distribution versus temperature. The general 1/y/p behaviour is probably due to 
mixing, The high temperature peaks reflect the initial condition, where the right 
cutoff moves quickly to lower temperatures again due to mixing. The peak around 
14,000 K is produced by a dynamical equilibrium between shock and compressional 
heating and cooling. It is located at a sharp rise in the cooling function. 





Fig. 4. Left: Energy loss in two Nirvana simulations with clouds and a control run. 
The control run as well as the 5c0 simulation are closed box simulations. Clearly, 
the cold clouds produce an enhanced system cooling rate. The clouds in the 6o0 
simulation are an order of magnitude denser, but it has an open upper boundary. 
Right: Cold gas mass fraction. The fraction of cold mass increases more rapidly 
with increasing cold mass content. This is not only true between the simulations 
but also for the time evolution of an individual one, which is an exponential growth 
till saturation at > 80% of the total mass. The control simulation produces no cold 
mass at all for all of the simulation time. 

100 Hence, if the velocity in jet cocoons is of order light speed, and the density 

101 ratio to the cooling b reak in Fig. 3 of order 10 4 , as suggested by other evidence 

102 (e.g. iKrauseL 2005bh . we derive velocities of the stirred up cool gas of order 

103 1000 km/s. This is indeed a figure that is often found for the emission line 

104 gas. 

105 We identify a prominent peak in the mass distributions over temperature at 
loe about 14,000 K (Fig. 3, right). This is not prescribed by the initial condition 
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107 as the other peaks, but due to a dynamical equilibrium between the turbulent 

108 energy transfer to the cool clouds and the cooling. Unsurprisingly, it is located 

109 at a sharp rise of the cooling function. 

no Energy loss and cold mass condensation increase drastically with increased 

in cold mass load (Fig. 4). The cold mass growth curve can be well fit by an 

n2 exponential. This is consistent with the naive expectation that the dropout 

n3 rate is proportional to the amount of cold gas present. 

H4 3 Discussion 



H5 3.1 Higher redshift radio galaxies 



lie As shown by Katherine Inskip at this meeting, emission line regions in con- 

H7 nection with higher redshift radio galaxies and the alignment effect are not 

us always within the region enclosed by the radio structures, which demonstrates 

n9 the complexity of the real situations. However, very often the gas is indeed 

120 within the radio structures. For such radio galaxies, the new mechanism pro- 

121 posed here may explain a number observations that were so far less well con- 

122 nected by physical theories. First, the observed emission line gas is steered 

123 up to high velocities, the order of magnitude of which can be explained by 

124 the Mach number - density relation. Second, the mechanism should produce 

125 a correlation between emission line and radio luminosity, since both is mainly 

126 powered by the turbulent e nergy in the jet cocoon. This is a well known 
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127 correlation in the literature (|McCarthvl . 11993). The mechanism may explain 



128 the increased visibility of the cold, optically-emitting gas within the radio 

129 structures just by the fact that the cold mass suffers a drastic increase due 

130 to the turbulence enhanced cooling. The observed evolution with source size 

131 (higher velocities and shock ionisation in smaller sources) may be understood 

132 if one considers that turbulenc e decays quickly. Turbulence in jet cocoons 
(FR Il lFanaroff and Rile v . 1974 ) is driven from the hotspots and decays away 



134 from the hotspots. As long as the source is small, the hotspots are closest to 

135 the galaxy where it can induce effectively multiphase turbulence. When the 
we hotspots are further away, the turbulence has already decayed when reaching 

137 the central parts where most of the emission line clouds are located. There- 

138 fore, the clouds are less stirred up, the velocity width is narrower, and there 

139 is less energy available for shock ionisation. Finally, the scenario also supports 

140 jet induced star formation. This is because as the cool clouds radiate their 
mi energy, they get more and more compressed and fragmented. The increased 

cloud density favours star formation. 



143 The simulations have to be confirmed in 3D which I am working on currently. 

144 Another issue is that the size of the effect may depend on the temperature of 
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145 the environmental gas. A change of the temperature of the gas surrounding the 

we host galaxies might then explain why the alignment effect appears at z > 0.5, 

147 only. 

us 3.2 Cold gas in nearby radio galaxies 



149 The jet-cold gas interactions reported bv lMorganti et al.l (|2005al ) are probably 

wo produced by a different situation than the one simulated here. However, since 

151 we end up with isotropic multi-phase turbulence, the model might be appli- 

152 cable in more general situations. The observations show ionised and neutral 

153 hydrogen co-spatial with the radio cocoon of 3C-305. The radio morphology 

154 suggests that the jet is hitting a disc and knocking off some cold gas. This 

155 seems to be accelerated to high outflow velocities, where both, neutral and 

156 ionised gas is observed. The simulations presented here also predict a close 

157 connection between neutral and ionised gas. In total, when the turbulence 

158 is evolved, there is a similar amount of T < 10 4 K gas present as gas with 

159 I0 4 K< T < 2 x I0 4 K, which might be taken as proxy for neutral and ionised 
wo hydrogen, respectively. The observations, however, suggest about a hundred 

161 times more neutral than ionised gas. The source of the discrepancy might eas- 

162 ily be with the simulations - a more realistic cooling function or 3D or maybe 

163 magnetic fields might make a difference. 

164 A connection to absorbers in high redshift radio galaxies, as suggested else- 

165 where, seems unlikely to me. Those are quite narrow whereas the one i n 3C 305 
is extr emely broad. The jet in wind model I have proposed earlier f Krausel . 



166 
167 



2005a[ ) seems not very appropriate here. The cold gas is probably not cooling 
lea behind a bow shock but entrained in a multiphase flow away from the radio 

169 hot spots. 
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